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' Abstract. We propose a new kind of supersymmetric (SUSY) transformation in the 

case of the two-channel scattering problem with equal thresholds, for partial waves of 
the same parity. This two-fold transformation is based on two imaginary factorization 
energies with opposite signs and with mutually conjugated factorization solutions. We 
call it an eigenphase preserving SUSY transformation as it relates two Hamiltonians, the 
scattering matrices of which have identical eigenphase shifts. In contrast to known phase- 
equivalent transformations, the mixing parameter is modified by the eigenphase preserving 
transformation. 
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^ ■ 1. Introduction 

The present work is a continuation of our previous investigations on supersymmetric 
(SUSY) transformations applied to coupled-channel problems with equal thresholds [1]. 
Our main aim here is to present a method based on SUSY transformations, which allows 
to construct potentials with given scattering properties, i.e., to solve an inverse scattering 
coupled-channel problem. 

There are several approaches to this problem based on the Gelfand-Levitan-Marchenko 
methods [2, 3]. In particular, Newton, Jost and Fulton [4]- [6] generalized the Gelfand- 
Levitan method and solved the corresponding integral equations in the case of two channels 
and rational scattering (S) matrices. Exactly solvable coupled-channel potentials obtained 
by this technique may be used for describing the neutron-proton scattering. In particular, in 
this way, Newton and Fulton [6] constructed a three-parameter phenomenological neutron- 
proton potential fitting low-energy scattering data. It would be interesting to extend this 
result by enlarging the number of parameters to fit scattering data on a wider energy range; 
however, the method based on integral transformations is rather involved and therefore 

t Boursier de l'ULB. 
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quite difficult to generalize. Using the Marchenko equation, the results of Newton and 
Fulton were nevertheless reproduced and improved by von Geramb et al [7]. A review of 
the present state of the art in the inverse scattering methods may be found in [8]- [11]. 

Our hope that the SUSY technique may be efficient for the multichannel Schrodin- 
ger equation is based on the well known equivalence between SUSY transformations and 
the integral transformations of the inverse scattering method for single-channel problems 
[12]- [15]. Due to this equivalence, one can use chains of first-order SUSY operators (also 
referred to as first-order Darboux differential operators [16]) for constructing a Hamiltonian 
with given scattering properties [17, 18]. This approach to the scattering inversion is more 
efficient [15] just because of the differential character of the transformation. There are 
several papers devoted to supersymmetric transformations for multichannel problems [19]- 
[26] (see also [27]- [30] for additional motivations and physical applications). Arbitrary 
chains of first-order SUSY transformations in the case of a matrix Schrodinger equation 
are studied in [31]. There, a matrix generalization of the well-known Crum-Krein formula 
is obtained. Another important ingredient of the supersymmetric inversion technique are 
the phase-equivalent SUSY transformations, which are based on two-fold, or second-order, 
differential operators. These are described in [32]- [34] for the single-channel case and in 
[24, 25] for the coupled-channel case. Such transformations keep the scattering matrix 
unchanged and simultaneously allow to reproduce given bound state properties. 

It should be noted that methods based on a direct generalization of the SUSY 
technique to the multichannel case are not able to provide an easy control of the scattering 
properties for all channels simultaneously. For instance, in the two-channel case, the S- 
matrix is parametrized by the eigenphase shifts 6±(k), 62(h) and mixing parameter e(k), 
where k is the wave number. Usual SUSY transformations modify these three quantities 
in a complicated way, which makes their individual control difficult. We believe that this 
is the reason why SUSY transformations did not find a wide application to multichannel 
scattering inversion. 

In the present paper, we propose a two-fold SUSY transformation which allows us to 
modify e(k) only, while keeping 5i(k) and 62(h) unchanged. We call such a transformation 
eigenphase preserving. It is necessary to stress the difference between this new kind of 
transformation and the well-known phase-equivalent transformations mentioned above. A 
phase-equivalent transformation does not modify the scattering matrix at all, whereas the 
eigenphase preserving transformation modifies the mixing between channels. An important 
consequence of that is the possibility to use single channel SUSY transformations to fit 
experimental values of the eigenphase shifts. Afterwards, the mixing parameter can be 
fitted without further modification of the eigenphase shifts. Thus, the main advantage 
of our approach consists in splitting the inversion problem into two independent parts: 

(i) fitting eigenphase shifts to experimental values independently for each channel and 

(ii) fitting the mixing parameter between these channels. To solve the first problem, one 
can use the single-channel tools mentioned above. In the present work, we propose an 
elegant solution to the second problem. 

In what follows, we will use definitions and notations introduced in our previous paper 
[1], where a first order coupling SUSY transformation is analyzed in details. Nevertheless, 
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in section 2, we recall some basic formulae necessary in the next sections. In section 3, 
we describe the new two-fold SUSY transformation and prove our main result that this 
transformation preserves the eigenphase shifts. A simple illustrative example of an exactly- 
solvable coupled-channel potential with a given scattering matrix is presented in section 4. 
In the conclusion, we discuss possible applications of the presented method and formulate 
some possible lines of future investigations. 



2. Two channel scattering with equal thresholds 

Consider the two component radial Schrodinger equation [35, 36] 

H Mk, r) = k 2 Mk, r), re [0, oo) (1) 
with the Hamiltonian 

#o = -1^ + Vj?\r) + 1(1 + IK 2 , I = diag(/ 1 , l 2 ) . (2) 

Here 1 is the 2x2 identity matrix and the interaction potential V mt (r) is a real and 
symmetric matrix, exponentially decreasing at large distances. We will consider the case 
of two partial waves li and Z 2 with identical parity, 

l 2 = k + 2m, m e Z. (3) 

For the sake of convenience, we combine interaction V^ int (r) and the centrifugal term into 
a single potential matrix Vo(r) = Vo nt (r) + 1(1 + l)r -2 . To characterize this potential near 
the origin, we use a matrix singularity index v. Matrix v is determined by the asymptotic 
behaviour of the potential near the origin, 

Vb(r^O) = v(v + l)r~ 2 + 0(1) . (4) 

Below, only potentials with singularity index being a diagonal matrix with integer entries 
v = diag(z/!, z/ 2 ) and Vj > lj are considered. We will call such potentials physical and 
restrict ourselves to SUSY transformations that produce physical potentials. 

As usual, the Jost solutions f (k,r) are defined as matrix solutions of (1) with 
exponential asymptotic behaviour at large distances [35, 36]. In what follows, we will 
need a more detailed asymptotic behaviour of these solutions; it is given by the asymptotic 
behaviour of the Bessel functions of the third kind, H\ l \(z), also called the first Hankel 
functions (see [37] for a definition). At large distances, the Jost solution thus behaves like 
the corresponding solution for the free particle 

i 

f (h,r^oo)^ding[h h (kr),h l2 (kr)} , h(z) = (^) 3 (z) (5) 

with 

fti(z-^oo) =6^1 + ^ + 0(^)1 , A = 1(1 + 1). (6) 
A special linear combination of the Jost solutions gives the regular solution 

Mk, r) = ^ [fo(~k, r)F (k) - f (k, r)F (-k)} , (7) 
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where matrix F (k) is the so-called Jost matrix. 

To construct eigenphase preserving transformations, we need solutions of the Schro- 
dinger equation (1) with a special behaviour both at large distances and near the origin. 
Thus, we first prove that the necessary solutions exist. 

Lemma 1. For any momentum k such that Im k > 0, det Fo(k) ^ 0, and for any constants 
Ci,2, d\,2 £ C, there exist two vector solutions u(k, r) and v(k, r) of the Schrddinger equation 
(1) which behave at large distances as 

u(k, r^oo)= e~ ikr ( Cl , c 2 ) T (l + o(l)) , (9) 

v(k, r^oo) = e ikr {d u d 2 ) T {l + o(l)) , (10) 

and near the origin as 

u(k, r ->■ 0) = (a ir Ul+1 , a 2 r U2+l ) T (l + o(r)) , (11) 

v(k,r ->■ 0) = (b 1 r--' /1 ,b 2 r~ U2 f(l + o(r)) , (12) 

where a^ 2 , 61,2 £ C 

Proof. To obtain the behaviour (10), r) can be expressed in terms of the Jost solution 

v(k,r) = f (k,r)(d 1 ,d 2 ) T . (13) 

Formula (12) follows from the behaviour of the Jost solution near the origin (see, e.g., [35]). 
Taking into account that Im k > 0, one gets from (7) 

ip (k,r -)• 00) -)• ^f (-k,r)F (k). (14) 

Here, we omit the second term in (7) since it becomes negligible at large distances with 
respect to the first term. Thus, solution u(k, r) may be obtained as 

Ik 

u(k,r) = —^{k,r)F^{k){ Cl ,c 2 ) T . (15) 
Formula (11) follows from (8). □ 

The 2x2 scattering matrix S (k) is expressed in terms of the Jost matrix as 

S (k) = e^F Q (-k)F -\k)e u ^. (16) 

Being unitary and symmetric, So(k) can be diagonalized by an energy dependent orthogonal 
matrix Ro(k) 

RT(k)S (k)R (k) = diag (e 2 ^ 1 ^), e 2iS °^ k) ) , (17) 

where 5o-j are the eigenphase shifts and the angle eo entering matrix R is called the mixing 
angle 

'co Se „W -sin £ „W 
' sme (A;) cose (A;) 
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Note that an opposite sign definition for the mixing angle could have been chosen; moreover, 
the order of the eigenphase shifts is arbitrary: exchanging them while adding ±7r/2 to the 
mixing angle keeps the scattering matrix unchanged. In the next section, the eigenphase 
preserving SUSY transformations are defined. 

3. Eigenphase preserving two-fold SUSY transformations 

3.1. Two-fold SUSY transformations 

Two-fold SUSY transformations lead to a number of interesting quantum models with 
unusual properties [38]. In particular, the corresponding superalgebra is nonlinear. It 
is natural to consider the two-fold SUSY transformation of the Schrodinger equation 
(1) as a chain of usual (i.e. one-fold) SUSY transformations. One-fold transformations 
for coupled-channel Schrodinger equations were introduced in [19]. Their generalization, 
which allows to introduce a coupling between channels, was given in [1]. The additional 
requirement that the transformed potential be physical was shown to result in a strong 
constraint on the transformation parameters. The case of two transformations is less 
restrictive since the intermediate Hamiltonian may be chosen unphysical. In particular, one 
may use as transformation functions complex-valued solutions of the Schrodinger equation 
corresponding to complex factorization constants. As we show below, a chain of two such 
transformations may preserve the eigenphase shifts. 

The chain of two SUSY transformations, Hq — > Hi — > H 2 , emerges from the following 
intertwining relations: 

L>iH = HiL\ , L 2 Hi = H 2 L 2 , (19) 

where the operators Lj map solutions of the Schrodinger equations to each other as 
■01 = LxipQ and ip 2 = L 2 ipi- These operators can be combined into an operator L defining 
the two-fold SUSY transformation 

LH = H 2 L, L — L 2 L 1 , (20) 

directly mapping solutions of the initial Schrodinger equation to solutions of the trans- 
formed Schrodinger equation as ip 2 = Lip . 

The operators Lj are first-order differential operators, 

L 1 = wi(r)-d r , L 2 = w 2 (r)-d r . (21) 

We use the standard notation for the superpotentials 

w j {r)=u' j {r)uj\r), j = 1,2, (22) 
w 2 {r) = u' 2 {r)u 2 \r), (23) 

which are expressed in terms of the matrix factorization solutions Uj and u 2 = L x u 2 . These 
solutions satisfy the following Schrodinger equations: 

H Uj = EjUj , E x u 2 = E 2 u 2 , (24) 
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with Ei, E 2 being factorization constants. Operator L then has a nontrivial kernel space, 
Ker L, spanned by the set of transformation functions Ui and u 2 . 

Ker L = span{-ui, M2} • (25) 
In the following, we will only consider self-conjugate factorization solutions, i.e. 
solutions with a vanishing self-Wronskian W[u, u] = 0. The Wronskian of two matrix 
functions u, v is defined as 

W[it, v\{r) = u T (r)v'(r) — u T \r)v(r), (26) 

leading for factorization solutions to 

W[ui,u 2 ](r) = u{(r) [w 2 {r) - wj(r)] u 2 {r). (27) 

Hence, self-conjugate solutions correspond to symmetric superpotentials. Solution u 2 then 
reads 

u 2 (r) = Liu 2 (r) = [w x (r) - w 2 (r)]u 2 (r) = - [v% (r)] 1 W[u ± , u 2 ](r) , (28) 

where the last expression has been obtained using (27) and the symmetry of w\. Using the 
Schrodinger equation twice, one also sees that the derivative of Wronskian (27) reads 

W[ Ul , u 2 ]'(r) = (E 1 - E 2 )ul{r)u 2 {r) , (29) 

a relation which will be used below. 

The Hamiltonians in (19) correspond to potentials related to each other through 
superpotentials 

V 1 (r) = V {r) - 2w' 1 (r) , V 2 (r) = Vi(r) - 2w' 2 (r) . (30) 

The sum of the two superpotentials W\ and w 2 defines the two-fold superpotential W 2 , 
which directly connects V to V 2 : 

W 2 {r) = Wl (r) + w 2 {r) , V 2 {r) = V (r) - 1W 2 {r) . (31) 

Using (23), (27) and (28), one can rewrite W 2 in the compact forms 

W 2 (r) = (E 1 - E 2 ) [w 2 (r) - ^(r)]" 1 (32) 

= {E 1 - E 2 )u 2 (r)W[ Ul , u 2 \-\r)ul{r). (33) 

As will be seen below, the second expression is more general than the first one, as it may 
be used in cases where the individual superpotentials w\ or w 2 are singular. 

Similarly, expressing the second derivative of the matrix solution ipo{k, r) from (1) and 
defining the logarithmic derivative 

w k (r)=^ (k,r)^\k,r), (34) 
one can rewrite the action of the second order transformation operator L on ip (k,r), 

ip 2 (k, r) = (w 2 - d r ) (wi - d r ) tp (k, r) , (35) 
in the following form 

Mk,r)= [{-k 2 + E l )l + W 2 {r){w l -w k )]^{k,r). (36) 
A more symmetric form of this formula 



ip 2 (k,r) 
may also be useful. 



,0 E 2 + E 1 \_ t TTr ,,fwi + w 2 
-k 2 + )1 + W 2 (r)[ w k 



M^r) (37) 
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3.2. Main theorem 

Let us now particularize the above results to two consecutive SUSY transformations with 
mutually conjugated complex matrix factorization solutions corresponding to imaginary 
factorization energies. We will prove that such a second order transformation modifies the 
mixing parameters without affecting the eigenphase shifts. 

Theorem 1. Consider a complex matrix solution u of the coupled- channel Schrddinger 
equation (l)-(4), with imaginary energy E\ — k\ = 2ix 2 and complex wave number 
ki = x{i + 1); X > 0; behaving at large distances as 

t v \ ( h h (-hr) ±ih h (hr) 
U{r °°' (-M MM I ' (38) 

and near the origin as 



«(r->0)=( * +1 * )[l + o(r)]. (39) 



a 2 r V2+l b 2 r~ 

The two-fold SUSY transformation defined by (20)-(24) with matrix factorization solutions 
Ui = u, u 2 — u* corresponding to the imaginary factorization constants E 1 , E 2 — E* — 
—2ix 2 and complex wave numbers k\, k 2 = x{i ~ 1); possesses the following properties: 

A. The resulting potential V 2 defined in (31) is real, symmetric and regular Mr. The two- 
fold superpotential W 2 reads 

W 2 (r) = 4ix 2 [w*(r) - w(r)]" 1 , w(r) = u'(r)u~ l (r), (40) 
= 4ixV(r)W[M,M*]~ 1 (r)M T (r), (41) 

where only the second expression can be used when the superpotential w is singular. 

B. The long range behaviour ofV 2 , 

V 2 (r ^oo) = 1(1+ l)r- 2 + o(r- 2 ), 1 = diag(Z 2 , h) , (42) 

corresponds to a re-ordering of partial waves with respect to channels. 

C. The scattering matrix S 2 of the transformed Schrddinger equation is expressed from the 
initial scattering matrix S as follows: 

S 2 (k) = O(k)S (k)O T (k) , (43) 

where the real orthogonal matrix O reads 

° (i) = e VMS - 2 * f H (44) 

D. The eigenphase shifts of the transformed scattering matrix S 2 coincide with the initial 
ones. With the permutation 

<M*0 = <M*0> (45) 
<M*0 = M*0> (46) 

the mixing parameter transforms as 

e 2 (k) = e (k) ± (-l) m arctan^-. (47) 

%X 
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Proof. First, we note that Lemma 1 implies that solution u exists. It reads 

«(r) = ^o(fci,r)F - 1 (fc 1 )^. °)+/o(*i,0^° * j . (48) 

Using (6) and (38), one may write the leading terms of the asymptotic behaviour of this 
factorization solution as 

( -ifcir ( 1 _ iAi \ ,1-Apikir ( i i iAi \ \ 
Tie-*- (l " ^) + j 

A. According to the choice of transformation functions and factorization constants, the 
one-fold superpotentials W\ and u>2 are mutually complex conjugated, W\ = w, W2 = w*. 
Therefore, one can use w = u'u~ l and its complex conjugated form w* in (22), (23) and 
(28), thus obtaining 

w 2 ( r ) = w*(r) = (u*)'(u*) -1 , u*(r) = LiU*(r) = (w - w*)u* . (50) 

In this case, (40) and (41) directly follow from (32) and (33). 

From (40), it is seen that W2, and thus the transformed potential (31), are real. The 
symmetry of matrix V2 (i.e. V 2 T = V2) follows from the symmetry of superpotential w, 
which can be established by considering the self-Wronskian W[iz, u\. Since (29) implies 
that this self-Wronskian is constant with respect to r and (49) implies that it vanishes at 
large distances, W[m,m](oo) = 0, one has W[w,w](r) = 0,Vr. According to (27), this is 
equivalent to the symmetry w T (r) = w(r),Vr. 

Let us now prove that V2 is regular. According to (31) and (41), this is the case if 
and only if the Wronskian W[u, u*] is invertible Vr. From (26) follows that W[u, u*] is an 
ant i- Hermit ian matrix, i.e. W[u, u*] = —W'[u,u*]. Moreover, using (29), the derivative of 
this Wronskian reads 

W[u, u*]'{r) = Aix 2 u T (r)u*(r) . (51) 

Its diagonal entries can thus be integrated using (39) and (49) respectively. One gets finally 



W[u,u*](r) = 



( Ai X 2 ]i\u ll {t)\ 2 + \u2i(t)\ 2 )dt W 12 [«,«*](r) 







-WI 2 [«,«1(r) -Ai X *j {\u l2 {t)\ 2 + \u 2 2{t)\ 2 )dt 

\ r / 



(52) 



where Uij and W^- [u, u*] label the entries of the factorization solution and of the Wronskian, 
respectively. This result implies that detW[w, u*] > 0, Vr, which proves the regularity of V 2 
stated in the theorem. Let us stress that this proof holds even in cases where superpotential 
w and the intermediate potential V\ are singular, which shows that expression (41), though 
more complictaed, is more general than (40). 

B. Let us first consider the case l\ 7^ / 2 . From the asymptotic behaviour (49), it 
follows that the determinant of the transformation solution u tends to zero as r — > 00 like 
the Laurent series 

det u(r -» 00) = (A ' - A ^ + o(r~ 2 ) . (53) 
X\l-i)r 
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Hence, the superpotential w behaves asymptotically as 

w(r^oo) = -^-( I ± M+0(1), (54) 




from which, using (40), we find the asymptotic behaviour of W 2 , 

W 2 (r oo) = _ i)+ o ( r_1 )- ( 55 ) 

It should be emphasized that from (55) follows the exchange of the centrifugal terms in V 2 
with respect to V [see (31)]. This effect of coupling SUSY transformations was previously 
described in [1]. Note that the scattering properties of the transformed system crucially 
depend on the exchange of centrifugal terms because of the presence of /-dependent factors 
in the S'-matrix definition (16). 

In the case of coinciding partial waves, l\ = l 2 , (55) is still valid but cannot 
be established through (54): instead, W 2 {r) can be calculated from the Wronskian 
representation (41) (see Appendix A). The fact that the two- fold superpotential vanishes 
at large distances faster than r _1 implies that the centrifugal tails are not affected by the 
SUSY transformations. 

C. To establish the modification of the scattering matrix, we have to look at the way 
the Jost solutions and the regular solutions transform in the two-fold transformation. 

Once again, let us start with the simpler case l\ ^ l 2 . Without loss of generality we 
may apply the general transformation of solutions (37) to the Jost solution, which now 
takes the form 



Lfo(k,r) = 



-ei + W 2 (r)[^^-w k 



f (k,r) = U(k,r)f (k,r). (56) 



2 

As we will see below, the matrix Uoo(k) = lim^oo U(k,r) determines the transformed 
Jost and scattering matrices. Using (54), (55) and the fact that W 2 Wk vanishes at large 
distances, one obtains a simple expression for this matrix, 



= I -e ) ■ < 57) 

From the dominant term of (5) and (6), it follows that the function 

f 2 (k, r ) = Lf (k,r)U^(k) (58) 

is the transformed Jost solution. 

As in the previous part, the case l\ = l 2 requires additional attention since the 
product W 2 (w + w*) gives at large distances the uncertainty • oo. Again we use the 
Wronskian representation (41) of the two- fold superpotential W 2 and the asymmetrical 
form of transformation (36) thus obtaining 

Lf (k,r) ->■ \(-k 2 + 2i X 2 )l + 4:ix 2 u*W[u,u*}- 1 u T '] f (k,r) . (59) 

Using (A.l) and (A. 6) in this expression leads to the same matrix U^k) as in (57). 
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Let us now find how the SUSY transformation modifies the behaviour of the potential 
at the origin. From (39), one gets 

det u(r -> 0) -> ai b 2 r ul ~ U2+1 - a 2 b 1 r V2 - yi+1 , (60) 

which suggests that the discussion will depend on the relative values of v\ and v 2 . 

For v 2 = v\, excluding the case a\b 2 = a 2 b\ (which requires higher order expansions), 
one can expand the superpotential w(r) in a Laurent series near r = 0, 

W (r ->()) = - [ aMUl + + a2blUl ~ a MM + 1) V q(1) (61) 

(aib 2 - a 2 bi)r \ a 2 b 2 (2u 2 + 1) -a 2 b x {v 2 + 1) - a x b 2 v 2 J 

which implies with (40) that the lowest-order term in W 2 is linear in r. Consequently, (31) 
implies that the singularity indices are not modified by the two-fold SUSY transformation. 
Note however that (30) implies that the intermediate potential V\ displays in general off- 
diagonal singular terms at the origin. 
For v 2 > vi, one gets instead of (61) 

w(r^0) = ±^ + 1 _° 2 j+o(l)- (62) 

To find the behaviour of W 2 at the origin, a higher-order expansion would thus be necessary. 
It is simpler in this case to study the two first-order transformations separately. From (30) 
and (62), we conclude that the intermediate potential V\ has the following singularity 
indices v — >■ v = diag(i/i + 1, v 2 — 1). For v 2 < u x , one gets v — > u = diag(^i — l,v 2 + 1) by 
symmetry. 

Let us now analyze the behaviour of the transformation function u* = L\u* which 
determines operator L 2 . Using (21) and (62) [or (61) when v x = v 2 \ one can find that a 
regular /singular vector solution transforms into a regular /singular vector solution of the 
new equation. Such transformations are called conservative SUSY transformations [26]. As 
a result the behaviour of u* near the origin is given by the conjugate of (39) with different 
values of constants a* l 2 and b\ 2 , i.e., a* l 2 — > a\ 2 and b* l 2 — > fo* 2 , and shifted singularity 
indices v = diag(z/i + 1, v 2 — 1) (to fix ideas, we consider the case v 2 > u x ) 

(~* 1/1+2 lf r -v x -\ \ 
kr- 4— 'J [1+0(r)1 ' <63) 

We have to split the discussion into two subcases, once again. For u 2 = z>i, i.e. 
v 2 = v\ + 2, an equation similar to (61) implies that w* behaves like r _1 multiplied by 
a non-diagonal matrix close to the origin. Consequently, the final potential V 2 will be 
unphysical in general, with non-diagonal singular terms at the origin; therefore, we will 
not consider this case any further. For v 2 > Vi, i.e. v 2 > V\ + 2, the same reasoning 
as above implies that the transformed potential V 2 has the following singularity indices: 
v — > v = diag(i>i + 1, v 2 — 1) = diag(z/i + 2, v 2 — 2). Finally, for v 2 < vi, which is the 
case for v 2 = v\ + 1, the second transformation restores the initial singularity indices 
v v = diag(i>i - 1, v 2 + 1) = diag(z/i, v 2 ). 
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The modification rules for the singularity indices of the potential may thus be 
summarized as follows in the physical cases: 

(1/1,1/1) -A (^1,^1) , (64) 

(1/1,1/1 + 1) A Ki/i + l), (65) 

(1/1,1/1 +m) -A (1/1 + 2,i/i + m- 2), m>2. (66) 

From here it is seen that in all cases Tr v = Tr £/. 

We are now ready to construct the regular solution of the transformed Schrodinger 
equation. For v 2 7^ v\ superpotentials w and w* have the structure given by (61) or (62) 
depending on the singularity indices. Therefore the first-order transformations Li and L 2 
are conservative. Thus, the result of the two-fold SUSY transformation applied to tpo(k,r) 
in the most general form can be written as follows 

L(p (k, r) = ip 2 (k, r)U (k) , (67) 

where Uq is a constant matrix with respect to r. Matrix Uo(k) is invertible V/c 7^ ki i2 , 
which can be seen from (25). In the case v 2 = v\, the conservativeness of the two-fold 
SUSY transformation can be established by considering (56) where ipQ is replaced by a 
regular solution. Note that <p 0t2 (k,r) = ipo t2 (—k,r); therefore, matrix Uq is an even matrix 
function of wave number k, U (k) = U (—k). The precise value of Uq is not important for 
the following. 

Applying operator L to the relation (7) between the Jost solutions and the regular 
solution, one obtains with (58) and (67) 

ip 2 (k,r)U (k) = ^ [f2(-k,r)U oo (-k)F (k) - f 2 (k,r)U oo (k)F (-k)} . (68) 
The transformed Jost matrix thus reads 

F 2 (k) = Uooi-tyFoWUoW . (69) 
The transformation of the scattering matrix then follows from its definition (16), 

S 2 (k) = e^U^e'^So^e-^U-'ik^ , (70) 

and is equivalent to (43) and (44). Note that the transformed S'-matrix does not depend 
on Uq. To prove that matrix O is real and orthogonal, one has to remember that li, l 2 , /1, 
l 2 all have the same parity, as implied by (3) and (42). If written like in definition (18), 
matrix O corresponds to a rotation angle 

2V 2 7T k 2 

T(-l) m arctan = - ± (-l) m arctan — . (71) 

D. Diagonalizing S2 in the same way as So in (17), 

R T 2 {k)S 2 {k)R 2 {k) = diag (e 2lS ^ k \ e 2lS ^ k) ) , (72) 

and taking into account that matrices Rq, O and R 2 = ORq all belong to 5*0(2), one sees 
that So and S 2 have the same eigenvalues. The mixing angle of S 2 is given by the sum 
of eo and (71). Inverting the order of these eigenvalues (see discussion following (18)), 
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one gets (45) and (47), i.e. a modification of mixing parameter vanishing at zero energy, 
e 2 (0) - e (0) = 0. □ 

3.3. Iteration 

Let us finally note that the transformed potential V 2 can be used as a starting point 
for a next eigenphase preserving transformation. This means that the two-fold SUSY 
transformation considered above can be iterated as long as desirable. A chain of n such 
transformations over the initial potential V will lead to the following mixing parameter: 

™ p 

e 2n (k) = e (k) ± (-l) m J2 arctan Yl (73) 

3=0 X i 

leaving the eigenphase shifts unchanged. 
4. Example 

Let us consider a simple example where the eigenphase preserving SUSY transformation 
is applied to an s — d diagonal potential with the following scattering matrix: 

n n\ i- A (k + iKi)(k + iK 2 )\ 

So(k) = d,a.g (l, L _ -A- _ .J j . (74) 

The first channel corresponds to the d wave and the second channel corresponds to the s 
wave. The corresponding potential reads 

6 

The rf-wave potential is purely centrifugal, while the s-wave potential is obtained from the 
zero potential by a second order one-channel SUSY transformation with the factorization 
solutions i>i(r) = sinh(/tir) and v 2 {r) = sinh^?")- This s-wave potential has no bound 
state but a singular repulsive core at the origin [17]. Potential Vo is thus characterized by 
the singularity and centrifugal indices 

v = diag(2, 2) , 1 = diag(2, 0) . (76) 

The Jost solution corresponding to potential Vq reads 



V (r) = diag^,-2(lnW[t; 1 ,t; 2 ](r)y / j . (75) 



where 



f (k,r) = diag(7o,KM,/o s (M), (77) 



/^)=e*(l + |-^), (78) 
/*(*, r) = ~ 9r V^t4 - 9 r ) e ikr N 1 N 2 , (79) 



v 2 {r) J \vxir-) 

with v 2 = [(lavi)' — d r ]v 2 and the normalization constants Nj = (ik — Kj) -1 . The regular 
solution <y9 is expressed from (7) with the Jost matrix 

F (k) =diag (1,-^3). (80) 
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2468 10 2 4 6 8 10 2468 10 

r > r 

Figure 1. Entries of the exactly solvable potential matrix Vi obtained from the uncoupled 
potential (75) with parameters Ki = 0.232 , k 2 = 0.944 by application of the eigenphase 
preserving transformation with x = 1-22. 



Using these expressions for the Jost and regular solutions, one may construct with (48) 
a transformation solution u with asymptotics (38) and (39), according to Lemma 1. The 
eigenphase preserving transformation described in Theorem 1 leads to a singular potential 
V 2 without bound state and with 

v = diag(2, 2) , [= diag(0, 2) . (81) 

The eigenphase shifts of the transformed S'-matrix coincide with the initial eigenphase 
shifts, 

k k 

S s (k) = — arctan arctan — , (82) 

Ki K 2 

5 d (k) = 0. (83) 

The mixing angle is given by (47) with eo = 0. In this case, different signs in (47) correspond 
to different signs in the coupling interaction V 2s d —V 2s( i- 

The transformed potential V 2 with the following parameters: 

ki = 0.232, k 2 = 0.944, x = L22 > (§ 4 ) 

is shown in Figure 1 [for definiteness we have chosen the positive sign in (47)]. The main 
reason to consider this example is that it illustrates the same scattering matrix as the 
one obtained by Newton and Fulton in [6]. The Newton- Fulton potential differs from the 
potential constructed above because it has one bound state. This difference can in principle 
be eliminated by the well known technique of the coupled- channel phase-equivalent bound 
state addition [25]. 

5. Conclusion 

In this paper, we have introduced an "eigenphase preserving" two-fold SUSY transfor- 
mation for the two-channel Schrodinger equation, i.e. a transformation that alters the 
mixing parameter between channels without modifying the eigenphase shifts. Chains of 
such transformations lead to coupling between channels in the scattering matrix which 
correspond to nontrivial fc-dependences of the mixing angle (73). With a reasonably small 
number of parameters, such mixing angles are probably able to fit experimental data, in 
a similar way to the usual phase shift fitting used in one-channel SUSY inversion [24, 18]. 
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Combining both techniques, we obtain a complete method of coupled-channel scattering 
data inversion based on SUSY transformations. As a first application of this method, we 
plan to invert the two-channel neutron-proton scattering data, hence improving the result 
of [6]. 

We also plan to study the following questions, raised by the present work. How do 
bound states transform under this eigenphase preserving transformation? How to construct 
a similar transformation for an arbitrary number of coupled channels? Do other forms 
of eigenphase preserving transformations exist? How will the presence of the Coulomb 
interaction modify the properties of the eigenphase preserving SUSY transformation? 
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Appendix A. 



Let us calculate asymptotics (55) using the Wronskian representation (41). This allows us 
to avoid manipulations with singular quantities which appear in (54) when l\ — l 2 . It is 
convenient to rewrite the asymptotic behaviour of the transformation solution in the form 



u(r -> 00) -> 2Q T - —KQ^a, 



Q T = (1 To- y )/2, £i = hr, 



(A.I) 



where A = diag(Ai, A 2 ), a x , a y and a z are the Pauli matrices, and the projection matrices 
Q T satisfy 

Ql = Q T , Q±Q T = 0, Ql = Q±, (A.2) 

Q±cr z = o- z Q T , Q±cr x = a x Q T . (A. 3) 

Here and in what follows we will only retain terms of order r _1 or lower. Let us first 
calculate the Wronskian asymptotics at large distances. Definition (26) leads to 

(1 



W(r ->• 00) ->■ Ai X e-^ iaz (a z ± a x ) 
which can be inverted (up to r _1 ) to give 



4 xr 



■(Ai + A 2 )a 2 Q q 



, (A.4) 



W-\r ->■ 00) ->■ — e-^ az 
8i X 



1 + 



4 xr 



-(Ai + A 2 )o~ z Q 1 



8i X 



1 



a z ±a x + — (A 1 +A 2 )Q ± 
2 X r 



(a z ± <7 x )e* lff *(A.5) 
^ . (A.6) 
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We can now calculate the two-fold superpotential up to r 1 

W 2 = 4ixVW -1 u T (A.7) 

-> x (^Q±a z (o- z ± <t x )Q± + — (Ai + A 2 )Q± - yQ±{<J z ± o- x )a z Q±Aj , (A.8) 
\>i "XT' si / 

where (A. 2) and (A. 3) have been used. To further simplify this expression, we also use the 

decomposition A = l(Ai + A 2 )/2 + a z {A\ — A 2 ), which leads finally to 

W 2 (r oo) ^(A 2 - A,)a z . (A.9) 

This expression is valid for any h and l 2 ; it is thus also valid for the case of coinciding 
partial waves. 
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